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FOREWORD 


This report sunnnarizes the major technical findings made in the 
research program at the University of Florida sponsored by NASA Lewis 
Research Center under Grant No. NSG-3018. This grant is a continua- 
tion of the support began under the same grant number on June 24, 3974. 

We report here our findings during the year ending September 1977, 

The work to be reported came from the efforts of several key 
people. C.T. Sah suggested to NASA the possibility of placing this 
grant research at the University of Florida, and afterwards, together 
with the author, helped make the arrangements and formulate the technical 
plan of the research. Dr, Sah has contributed significantly to nearly 
every aspect of the research described in this report. After joining the 
faculty of the University of Florida in September, 1976, Arnost Neu- 
groschel became a key contributor to the research program. The technical 
collaboration of M.P. Godlewski and H.W. Brandhorst, Jr., of NASA Lewis 
Research Center greatly aided the completion of several phases of the 
research. During the tenure of the NASA grant, the Energy Research ?nd 
Development Administration began supporting a separate but related re- 
search program at the University of Florida under the author’s direction. 

This support aided the research reported here, and made possible a 
valuable and continuing technical interaction with J.G. Fossum of Sandia 
Laboratories (ERDA), Albuquerque, New Mexico. In addition to these senior 
investigators, two graduate students at the University of Florida participated 
in the research: P.J. Chen and S.C. Pao. 
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CHAPTER 1 INTRODUCTION 


As Brandhorst first noted, the power conversion efficiency of 
p-snbstrate, p-n junction silicon solar cells is limited because Vqq 
falls short of theoretical expectations by about 15 percent [1]. Thus, 
increasing V „ is important in raising efficiency, doubly so because 
the fill factor FF increases with increasing [2]. Indeed, as we 
have recently pointed out [3], essentially in agreement with the earlier 
projection of Brandhorst [1], if were increased to 700 mV, AMO 
efficiencies in the range of 18 to 20 percent become possible. 

1.1 Studies of Basic Mechanisms 

Our previous work at Florida under NASA LeRC grant support has 
proposed methods, based on experiment [4], for determining which of the 
many fundamental mechanisms [5,6] that can occur in silicon devices are 
responsible for this discrepancy in in any given solar cell. This 
work has demonstrated experimentally the important conclusion that the 
discrepancy originates in the excessive recombination current occurring 
in the degenerately doped emitter of the cell. In particular, for cells 
having a base resistivity of about 0.1 Ohm-cm, we demonstrated that the 
dark recombination current from the quasineutral emitter region exceeded 
that from the quasineutral base region by more than an order of magnitude 
[4], We demonstrated further that both of these recombination currents 
exceeded the dark recombination current coming from the junction space- 
charge region. In these demonstrations, the dark cell was biased in the 
forward direction with a voltage V near 600 mV, which is approximately 
the value of V^^ produced in these cells. The conclusion implied by 
these experiments is that the large recombination current occurring 
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either in the bulk of the degenerately doped emitter or at its surface 
limits V and efficiency well below the projections of classical p-n 
junction theory [1]. This conclusion, based on experiment, contrasts 
with the recent theoretical suggestions of other authors that field- 
induced bandgap narrowing in the junction space-charge region is the 
dominant contributor to the observed performance, including and 
efficiency, of p-substrate, diffused p-n junction silicon solar cells 
[7,8]. 

As yet, the physical mechanisms causing the anomalously large 
emitter recombination current have not been identified. However, 
considerable progress has recently been made. Our work at Florida has 
shown experimentally [9] the decided superiority of forward-biased 
admittance measurements over open-circuit-voltage-decay (OCVD) measure- 
ments in detertaining the dominant relaxation time [4], originally called 
the dominant natural mode [10], of the cell. The accurate determination 
of the dominant relaxation time is crucial to the determination of the 
minority— carrier charge stored in the quasi-neutral emitter and of the 
phenomenological carrier lifetime of the emitter, and thus to the 
determination of the physical mechanisms responsible for the excess 
emitter recombination current [4]. Furthermore, we have just reported 
a new method, based on the temperature dependence of the emitter recombina- 
tion current [11], for determining the bandgap narrowing in the quasi- 
neutral emitter region as a function of the emitter doping concentration. 
Bandgap narrowing is believed to be one of the main physical mechanisms that 
contribute to the discrepancy in [4,5,6]. This new experimental method 
developed at Florida, together with work just described by Lanyor. and co- 
workers 112] and by Mertens and co-workers [13], may yield insight needed 


3- 


to better define the fundamental physical mechanisms underlying the ^OC 
discrepancy. 

1.2 Studies of HLE Structures for Increasing 

In the past year, although our work on fundamental physical mechanisms 
has continued to some degree, the main attention of the NASA LeRC sponsored 
research at Florida has shifted toward the demonstration that new device 
structures can substantially diminish the effects of whatever mechanxsms 
in the emitter may be degrading without requiring a detailed 

Identification of these mechanisms. This new emphasis led to our proposal, 
from a theoretical standpoint [14,15], of the high-low-junction emitter 
(HLE) solar cell, in which the emitter region contains a high-low (H-L) 
junction. We have demonstrated experimentally [3,16] that this H-L 
junction can suppress the dark emitter recombination current to levels 
so low that the emitter current no longer presents an obstable to the 
achievement of an open-circuit voltage of 700 mV. For test cells we have 
fabricated at Florida, the open-circuit voltage observed has been mainly 
determined by the dark base recombination current. Our experiments have 
shown that large values of V result in these HLE test cells, but that 
values approaching the theoretical limit of 700 mV x^ill require improving 
the diffusion length of electrons in the quasineutral base region. 

The maximum V _ seen to the present has been 640 mV for AMO illumina- 

oc 

tion. This voltage, which compares with the 600 mV maximum normally seen 
in p-n junction silicon solar cells, has been measured on an oxide-charge- 
induced (OCI) HLE structure. In this structure, the H-L junction is 
formed in emitter material of non-degenerate doping concentration by an 
oxide-charge-induced electron accumulation layer [3,16]. Open-circuit 
voltages about 10 mV less have been observed in another type of HLE solar 
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cell, the diffused-HLE structure, in which impurity diffusion forms the 
H-L junction in the emitter region. 

1,3 Structure of This Report 

This report consists of three parts. In the first part, which Is 
Chapter 2, we develop, from a theoretical standpoint, the modifications 
of the basic Shockley equations that result from the random and non- 
random spatial variation of the chemical composition of a semiconductor. 
These modifications underlie the existence of the excessive emitter 


recombination current that limits the of solar cells. The discussion 
includes a more rigorous treatment of quasi-fields than has appeared 
before. In the second part, work related to the measurement of parameters 
is presented. Chapter 3 treats the measurement of series resistance. 
Chapter 4 deals with the measurement of the base diffusion length and 
presents two methods for establishing the energy bandgap narrowing in 
the heavily-doped emitter region. The discussion in Chapter 5 relates to 
corrections that can be important in the application of one of these 
methods to small test cells. In the third part, we briefly describe 
OCI-HLE test cells which exhibit a considerably higher ^OC than was 
previously seen in n-on-p solar cells. 


Insofar as is possible, each of the Chapters are written as in- 
dependent uiiits to enable the reader with a main interest in a later chap- 
ter to be able to read that one directly without first having to become 


familiar with the earlier chapters- 
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CHAP. 2 CARRIER GENERATION, RECOMBINA'I ION, TRAPPING AND TRANSPORT 
IN SEMICONDUCTORS WITH POSITION DEPENDENT COMPOSITION 


2.1 Introduction 

Semiconductor structures and materials with position dependent chemical 
composition have been of central interest in device applications. For structures 
with abrupt changes in but otherwise uniform chemical compositions, such as 
the hetroj unction Interfaces of a Schottky barrier or a MOS capacitance, the 
mathematical analyses of carrier generation, recombination, trapping and 
transport (GRTT) are based on ideal and position independent energy band 
models for each material region [1,21. The solutions of the device charac- 
teristics are then obtained by matching the solutions of each of these uniform 
regions at their interfaces. These solutions are obtained from the following 
Shockley equations which were first used to get the characteristics of p-n 
junction diodes and transistors [3,4]. These were generalized to include 
the most important generation, recombination, trapping and tunneling mechanisms 
[5] and high-frequency effect [6] near thermal equilibrium as well as hot 
carrier effects on transport [7]. 
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0 = - VJp + 


7 . I = - = p/e = (q/E)(P - N + N 


DD “ *^M “ 


( 4 ) 

( 5 ) 


Here, the coefficients which characterize tiie material properties aiul the 
variables have the usual meaning [5J, 

When the chemical composition of the material or device structure varies 
continuously and strongly in a regular or random fashion, the Shockley equations 
need to be modified. A number of fairly distinct situations can be treated 
and discussed below. 


(1) If the composition varies on an atomic scale but the atomic perturbations 

17 3 

are widely separated (concentration less than about 10 perturbations/cra ), 
we have the classical case which is contained in the Shockley equations; the 
mobilities, p and p , the dif fusivities , D and D , and the generatlon- 
recombination-trapping rates all vary with position. These material parameters 
are macroscopic averages taken in regions of sufficiently large volumes to 
contain many (- 10 ^) particles (electrons, holes and atoms) to make the statistical 
fluctuation small. 

(2) If the composition variations are on a macroscopic scale (concentrations 

17 —3 

of the foreign atomic species or physical defects are greater than about 10 cm 
but smaller than about 1% of the host atomic density), the Shockley equations are 
still applicable. However, the ideal energy band model of allowed band of 
nonlocalized energy (Bloch) states or levels with sharp band edges separated by 
forbidden energy gaps must be modified. The modification comes about from 
the high concentrations of the spatially random atomic perturbations which shift 
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the band states of the ideal model Intu the energy gaps, thereby producing a 
distribution of semi-localized energy gap states both in energy and in position. 

The deep gap states, which are located far away (many kT/q in energy) from 
the ideal band edge positions, are nearly true bound states with negligible 
carrier mobility. They can be accurately included in the generation-recombination- 
trapping terms of the Shockley equations. However, the shallow gap states, 
which are located in energy near the ideal band edges, may have appreciable 
mobility due to the large spatial overlaps of the very extended wave functions. 
These edge states will modify the macroscopic and thermal averages of the 
mobilities and dif fusivities in the Shockley equations, giving them position, 
temperature and electric field dependences significantly different from those 
of the ideal energy band model discussed in case (1). Here, the ideal or 
sharp band edges lose their significance. To put the physics of this situation 
into a simple persepctive, Cohen [8] introduced the concept of mobility edge 
to separate the semi-delocalized edge states from the nearly localized tail 
states. The modifications of the macroscopic and thermal averages of the 
mobility and dif fusivlty from tlie ideal case, then, come in two factors; Che 
integration limit is now the position of the mobility edge instead of the ideal 
band edge and the density of state in the integrand is decreased due to the 
splitting of the band states i> o the gap. Van Overstraeten [9] has provided 
an analysis of Che mobility in this situation. Although he included the Call 
states in the mobility and diffusivity averages by assuming that tliey have the 
siimo carrier scattering power or tolaxnCion time as the band states, the errors 
in computing the diffusion and drift current are probably small since the densities 
of the electrons or lioios trapped at these localized gap states are much smaller 
Chan those in the edge and band states. 
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The effects of the distributed gni) states on the carrier generation- 
recombination rates are well known. It was demonstrated as tlie principal 
source of excess current in tunnel diodes [10]. However, a demonstration of 
the application of this model for calculating the recombination rates is 
not available. An example will be given in the next sectiont, 

(3) If the composition variations are on a chemical scale and continuous 
over an extended region, basic modifications of jhe Shockley equations must be made. 
These modifications can be analyzed by separating out the random and the 
nonrandom components of the atomic potentials. The random component will 
produce the modifications to the Shockley equations similar to those discussed 
in cases (1) and (2). The nonrandom and continuously varying component will 
produce a corresponding spatial variation of energy levels of the band or Bloch 
states as well as the semi-localized gap (edge and tail) states. This produces 
major modifications of the Shockley equation, giving rise to new terms. The 
position dependent energy level was first analyzed by Kroemer [II] for device 
applications in graded or position-dependent band gap materials. He began with 
the effective mass equation and introduced the concept of quasi-electrical 
fields for electrons and holes. 

Although the graded band gap materials have seen little application in 
conventional transistor devices in the past, it shows promise for improved 
efficiency in solar cells. Thus, an extension of the Shockley equation for this 
case will also be presented. 

In the next two sections, we shall obtain quantitative results for the 
giMUMMl ioii-rt'i'ombfmil Ion ratos in lUo dlstrlhntod gup slat es of case (2) and 
modifications of the Shockley equations in the graded band gap materials of 


case (3). 
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2.2 Generacion-Ilecombinatlon-Trapplri)', \n Distributed Gap States 

For this case, the basic property of the gap states is that each imperfection 
center, whether having a point, pair or more complex geometry, is separated from 
its neighboring centers by a distance large compared with extent of the bound 
state wave function localized at the center. Thus, the transition rate of an 
electron (or hole) bound to one center to the adjacent centers is much smaller 
than the rate to the conduction or valence band and edge states. Then, the capture 
(or emission) rates of electrons (or holes) at each gap-state center can be added 
to give the total rate. 

The above property of independent or noninteracting gap-state centers make 
the arithematics of the original Shockley-Read derivation [12] of the recombination 
kinetics directly applicable. If the centers were interacting, then additional 
rate equations for the bound- to-bound transitions between the centers or between 
states on one center must be Included, This was done previously [31. 

For the independent centers, the four electron and hole capture and emission 
rates between a band state and trap state located at are shown in Figure 1. 

The transition rates, averaged over the distribution in the conduction or valence 
band states are just those given by Shockley and Read: 

^T ^ (6) 

= e^(E,j.)f^,(E^,)P^(E^)dE^ , (7) 

d^ = ep(E^)[l-f^(E^)]P^OydE^, (9) 

In Equation (6) , is the macroscopic averaged capture probability rate 
(///cm ^-Hcc) of electrons in the conduction band by the gap state located at 
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an energy of In Equation (7), is the corresponding emission rate 

of electrons, n is the concentration of electron in the conduction band, 
f^(E^) is the fraction of the gap states occupied by electrons at ttie energy 
and p(E^)dE^ is the concentration of the gap states located in the energy 
range and E^+dE^- The coefficients In Equations (S)and (9) has the 
corresponding meanings for holes. 

Although the rate equations (6) to (9) are used by Shockley and Read 

specifically for the thermal transition process, they are also applicable to 

optical and Auger-impact transition processes using extended definitions for 

the capture and emission rate coefficients, c , c , e and e C5l. 

’ n’ p’ n p 

The total capture and emission rates through all the gap-state centers 
can be obtained by sunimiiag or integrating Equations (7) to (9) over the 
distribution of the gap states in the band gap. However, since the gap states 
are independent, we may apply the condition of steady state to Equations (6) 
to (9) before making the summation over the gap states. This will be used 
as an illustration of the modifications of the recombination rates by the 
presence of a distribution of gap states. Thus, under the d.c. steady-state 
condition, we have 

3t - - Ct + = (a^f^)/dt = 0 (10) 

which gives the well known Shockley-Read stp'dy-state results 




c^(E„)N + e^(E ) 
T JL ^ 




( 11 ) 
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( 12 ) 
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The total net steady-state recomb inat inn rate is then the sum of R^(E^) 
over all the gap states: 


K, 


^SS ^ 


c (Ejc_<EjKP - e_(Eje_(Ej 


r T 


r p' T' 


I c (E„)N + e (EJ + c (E^)P + e (E^) 
Ey n T n T p T p T 


p^(E^)dE^ 


( 13 ) 


Here, the integration is from the mobility edges, E to E_, over the density 

V u 

distribution of the gap states, p^(E^). 

Let us now work out a simple example for solar cell junctions which has 
an explicit result. Consider the forward bias condition in a p-n junction. 

The constant quasi-Fermi level approximation may be used in the junction space 
charge layer so that we have 

PN = n^exp(qV/kT) 

and if we neglect the spatial dependences, then 
P = N ~ n^exp(qV/2kT) 


(lA) 


(15) 


The steady-state recombination rate given by Equation (13) then simplifies to 


n^(exp(qV/kT) 


1 ) 


SS 


Cexp(qV/2kT) + 1 






(16) 


where the ratio C is somewhat dependent on the properties of the gap states as 
a function of energy but is assumed to be roughly unity in this example. It 
is unity if (c +c )n, = e +e at all gap-state energy levels. *.s an 

interesting simple result which shows that the junction current due to recombination 
in the space charge layer for a distribution of gap states is proportional to the 
geometrical average of the emission rates over the distribution of the gap states. 
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2.3 Position Dependent Chemical Composition 

For this case, we will consider only the effects from the nonrandom 
component of the atomic potentials since the gap states produced from the random 
part can be taken into account by the way described in section II. The 
fundamental starting point of this problem Is the many*'body Schroedinger 
equation for the entire crystal which is assumed to be large so that surface 
effects are unimportant. Due to position denpendence of the chemical 
composition, there is a slowly varying aperiodic crystal potential which is 
the sum of the potentials of all the atoms in the crystal. Then, in che usual 
manner, the core electrons attached to each atomic nucleus can be taken into 
account by a transformation [i3j to give an effective aperiodic pseudopotential. 
The aperiodic part can then be isolated out from the periodic part and treated 
as an extended perturbation to Che periodic potential problem. The pseudo- 
effective mass equation L14J can then be derived for the perturbed problem. 
Although a detailed mathematical analysis following the above procedure has not 
been worked out for the specific problem on hand, a simple physical picture 
emerges x^hich is certainly valid for sloxjly varying chemical composition. A 
quantatlve condition of validity is that the perturbation potential in an 
extended region under consideration should be small compared with the 
periodic electric field due to atomic Ions in that region. An equivalent 
condition is that the perturbation potential referenced to an average potential 
in the extended region should be small compared with the energy gap of interest. 
For such a slowly varying perturbation, the spatial varation can be represented 
by tlie position dependences of tlie electron affinity, X(r), and energy gap, 
x'r a I terimt IveJy by the position dependences of the conduction and 
valence band edges, and (r), which serves as the effective potential 

energies seen by the electrons and holes In the effective masH equal Iona. In 
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addition, there are also spatial depedences of the density of state or effective 
nt/innoH which muni h<‘ fnkon tiilo ticcottni , Tliln hand inculcl la ilhiHtrated in 
Figure 2. 

The equilibrium one-electron distribution function is then given by 

f(E) ^ (17) 

1 + exp(E-E^)/kgT 

where the total energy is the sum of the potential energy and the kinetic energy 
and E is the Fermi energy. In the conduction band, E > E and 

i? . Cj 

E = + (nV/2m^) (18) 

and in the valence band, E < E^ and 

E = (19) 

For nonequilibrium, we can define a quasi-Fermi level for electrons, F^, 

and one for holes, Fp, so that the distribution functions have the form of 

f(E) (20) 

1 + exp(E-E)/k„T 


where for electrons and F=Fp for holes. Using the standard perturbation 

procedure to solve the Boltzmann equation C6] , the modified current equations 
corresponding to the ideal case given by Equations (1) and (2) are 


J„ = + qD VN - qy NVV,, 
N n n C 

(21) 

= - qy NVV„ 
n N 

(2ia) 

and 



(22) 

= - qy^PVVp 

(22a) 

where V^,=-E^/q, Vj^=-Fj^/q, Vy=-Ey/q and Vp=-Fp/q. 

For the nondegenerate case, 


the electron and hole concentrations are given by 
N = N expC(E -F„)/k_T] 

C C N B j23) 

P - N,,expL(Fp-B^)/kjT] 

The mobilities and diffusivlties appearing in Equations (21) to (22a) are 
now the local values which are averaged over the carrier concentrations N and P. 
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In addition, the effective density oE states, and N^, are also local va Liioh 
which are spatially varying due to the position dependences of the density of 
energy states. 

The conventional concepts of an intrinsic sample as well as the intrinsic 

carrier density and intrinsic Fermi level at thermal equilibrium also lose the 

familiar meaning since they are functions of the spatial dependences of E (r) , 

E (r), E„(r) and example, the position dependences of the electrostatic 

potential for a pure semiconductor, i.e. a graded band gap semiconductor whithout 

any ionized chemical impurities and defects, can be obtained only by solving the 

Poisson equation in a given material with known spatial dependences of E^, 

E , E and X specified boundary conditions. This electrostatic 

potential variation with position, however, will not be the same as the spatial 

dependences of the intrinsic Fermi potential, E^{r), which can still be calculated 

in the conventional way using Equation (23): N=P=n.=/N~N exp[ (E„-E )/k T] and 

1 G V G V D 




The example just given seems unfamiliar at first glance since it defies our 
firmly-rooted, concepts used in semiconductor device analysis. However, 
further thought on a familiar situation immediately removes this conceptual 
bottleneck. For example, the d if ferent spatial variations between the intrinsic 
Fermi potential and the electrostatic potential in a graded band gap material 
has an analogous counterpart in the contact region between two dissimilar but 
compositional uniform materials such as the oxide-silicon interface in a MOS 
capacitor. In this case, if we let the electrostatic potential in the oxide 
coincide with the intrinsic Fermi potential of the oxide and disallow a 
tU Hconl, I lui i ly in the electrostatic potential at the Interface plane betw ,*'n 
the oxide and tlie semiconductor, then the electrostatic potential In the 
semiconductor will not coincide with the semiconductor intrinsic Fermi potent la J 
but t^ill be displaced from it. From this interface picture, one can then readily 
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build up the spatially varying band pirture shown in Figure 2 for a graded gap 
semiconductor by placing many macroscupLcally thin semicondcutor slices 
with different chemical compositions In contact. The conclusion on the 
difference between the electrostatic and intrinsic Fermi potentials is then 
immediately obvious. Such a multi-layer structure is in fact the model 
one implicitly uses to solve the one-electron Schroedinger equation as well as 
the transport equation In a one-dimensional graded gap semiconductor. 

A numerical demonstration of the differences between electrostatic and 
intrinsic Fermi potentials can also be inferred from the work of Emtage which 
was undertaken many years ago Cl5J. 

An important point to be noted is that the electrostatic field, given by 
the gradient of the electrostatic potential, loses its fundamental microscopic 
significance in a graded band gap semiconductor. It is the gradient of the 
conduction or valence band edges, VE (r) and VE (r) , which appears in the 

V# V 

Schroedinger equation and which give the forces experienced by the electron and 
hole quasi-particles and not the gradient of the electrostatic potential. This 
point was made clear by Kroemer Cll] who coined these effective fields the quasi- 
electric fields [16]. Thus, the hot electron properties in a high electric 
field is dependent on the quasi-electric field for electrons, VE^(r) , and not 
the ordinary electrostatic field, VV^., obtained from the Poisson equation. On 
the other hand, the boundary conditions at an interface between two dissimilar 
materials are to be imposed on the electrostatic field, electrostatic 
displacement and electrostatic potential or the true macroscopic charge dnesities, 
N, P, Njjp and but not on the quasi-electric fields. 
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2.4 Sununary 

The modifications of the Shockley equations by the position dependences of 
the atomic composition of semiconductors are described. Two modifications are 
analyzed. The first concerns the random spatial variations which give rise to 
distributions of highly localized states in the band gap. These band gap states 
enhance the carrier recombination rates. The total steady-state recombination 
rate at these gap states is worked out as an example and applied to the space 
charge region of a pn junction. 

The second modification discussed in this paper is a major one. It is 
shown that change in the form appears in the two current equations of the five 
Shockley equations for a two-band semiconductor. Instead of a proportionality 
to the conventional electric field, the electron and hole drift currents are 
now proportional to the gradient of the conduction and valence band edges 
respectively. These gradients, of band edges are also known as the quasi-electric 
fields and they are not equal in the graded band gap materials while they are 
equal in an uniform material. The quantum mechanical basis for the presence 
of the quasi-electric fields are discussed and shown to come from the aperiodic 

t 

part of the potential energy of the crystal in the many-body problem and its 
one-electron representation. The use of the quasi-electric fields for hot 
carrier effects and the electrostatic field for boundary condition matching 


are also pointed out. 
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Figure Captions 

Figure 1 The four electron and hole capture and emission transition 
processes between a conduction or valence band state and a 
bound state located at an energy 

Figure 2 

(a) The one-electron energy band model of a solid plate with position 
dependent composition giving rise to spatial variations of 

the energy gap, E_, electron affinity, x» the band edges of 

the conduction and valence bands, E and E . The physical 

V 

crystal potential composed of the sura of the atomic potentials 
is not shown in the figure. The shaded regions are the 
conduction and valence bands. Gap states such as those 
indicated in Figure 1 are not shown here for simplicity. 

The, rise of the potential well (heavy curves) near the crystal 
surface is highly expanded in relation to the width of the 
crystal labeled. 

(b) The one-electron energy band picture of two semiconductor crystals 

in contact giving rise to spatial variations of the energy gap, E (y) , 

electron affinity, X(y)> the conduction and valence band edges, 

Ep(y) and E (y) . The two vacuum levels, on the left and right, are 

shown to illustrate a large area semiconductor plate situation. The 

dashed E^Cy) and E^(y) curves with discontinuities are those of the 

idealized and physically unrealistic situation assumed in the past 

which overlooked the distortion of the atomic structure and electron 

distribution in the interface region which give the gradual transition 

of X» E_ and E„ shown in solid curves. 
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CHAPTER 3 EXPERIMENTAL DETERMINATION OF SERIES RESISTANCE 
OF P-N JUNCTION DIODES AND SOLAR CELLS 


3 . 1 Introduction 

The series resistance Rg is a parasitic element in diodes 
and solar cells that can degrade the device performance. In 
general/ the presence of R_ will bend the current-voltage 
characteristic away from the idealized characteristic [1] - 
In a solar cell/ this bending decreases the power conversion 
efficiency [2] . 

The series resistance of a diode or a solar cell consists 
of the contact resistances of the metal contacts at both faces 
of the device and of the resistance of the bulk semiconductor. 
The resistance of the bulk substrate can be calculated taking 
into account a spreading resistance correction [3]. The con- 
tact resistance is difficult to calculate accurately. In 
solar ceils with a very thin diffused layer and contact grid 
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on the top, the jresistance of the diffused layer might 
dominate. Calculation of this resistance is difficult 
because of the nonuniform current flow occurring in the 
diffused layer [4]. 

Thus, because of the uncertainties and difficulties 

involved in calculating methods for direct measurement 

of R are desirable. For a solar cell, a method has been 
s 

described [5] that involves exposing the solar cell to 
different illumination levels. Here we describe three 
other methods for measurement of R^, which can be applied 
to any diode or solar cell. These methods involve: 

a) small signal admittance measurement at 
moderately high frequency; 

b) combined measurement of dc and ac conductance; 

c) open-circuit-voltage -decay measurement (OCVD) . 

3 . 2 Theory 

The small signal equivalent model of the ideal junction 
diode, derived by Shockley [1], consists of a diode conductance 
and a diode capacitance in parallel. The circuit 

elements, and result from solving the continuity 

equation for holes and electrons under small signal excita- 
tions including both the quasineutral base and emitter re- 
gions. In parallel with these two elements, the circuit 
representation of a real diode includes two capacitances 
originating in the junction space-charge region; C^, 
associated with the ionized impurities, and 


associated 
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with the mobile charges {6], The total capacitance is 
Cp = Cgjq + + ^'SCR* complete the circuit, Rg is 

added, as shown in Figure 1(a). 

Method (a) i As shown by Sah [6] , increases exponen- 
tially with increasing bias voltage whereas, for small 
forward bias, - 0,^, shows a relatively weak dependence 
on Vq. If the small“signal admittance of the diode is 
measured at a forward bias small enough and at a frequency 
to high enough that o)Cj^>>Gj^, then the diode equivalent cir- 
cuit reduces to R^ and in series, as shown in Fig. 1(b) 

Rg can then be easily measured using a suitable high- 
frequency (1 to 10 MHz) bridge. 

Method (b) i In contrast to method (a) , we now set the for- 
ward bias voltage V and the frequency to to such values that 
toCp-Gj^. Noting that an admittance bridge measures Gp and 
Cp of the parallel combination shown in Fig. 1(c), we see that 
the dependence of Gp plotted against V will differ from 
the dependence of Gj^ against V because of the presence of 

R_. For dc conditions, as illustrated in Fig. 1(d), the 
s 

measured conductance Gj^^^ = [Rg + 1/^D^ ^ will show yet another 
dependence on V. Comparison of Gp(V,Rg) versus Gj^^(V,Rg) 
then yields R_. A simple computer program can be written to 
enable rapid determination of R^. 

Method (c) : The open-circuit-voltage-decay (OCVD) measurement 

[7] has been often used to determine minority carrier life- 
time of a diode. In this measuiement , a voltage step applied 
in the forward direction causes a flow of forward current I. 
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The terminals are then open-circuited. The resulting 
cessation of I causes a sharp drop of the voltage AV to 
occur at the diode terminals. This voltage drop AV 
divided by I gives series resistance R^. 

3 . 3 Results 

Table I shows results from methods (a) and (c) together 
w.ith data obtained by the illumination method of Wolf and 
Raushenbach [5], The devices used in the measurements were p -n 
diodes, 30 mil. in diameter and 2x2 cm n -p solar cell. The 
OCVD method was inaccurate for the devices measured since 
the voltage drop AV was difficult to determine accurately; 
this is especially difficult in devices with short life- 
times. Method <b) could not be used for our devices, since 
the ac measurements have to be taken at a frequency co such 
that 0) << 2/t [1], where x is the minority carrier lifetime 
in the base. In the devices we studied, x was a few micro- 
seconds which requires measurement at frequencies not ex- 
ceeding about 1 kHz. At such low frequencies, Gp will dif- 
fer from by only about 0.005%, which is much less than the 
estimated experimental accuracy of 1%. This method can be 
used, however, for short- lifetime devices. For example, 
if X = O.lysec, the ac measurements could be done at a fre- 
quency of a few hundred kHz, at which Gp - G^^^ would be 
large enough to enable a determination of R^. 

Method (a) is the simplest approach and gives accurate 
results that agree reasonably well with those produced by 
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the illuinination method. The accuracy of method (a) is 

limited by the accuracy of the bridges used, which is 

better than five percent for the measurements reported here. 

Measurements were taken at frequencies between 1 and 5 MHz. 

The results obtained from the illumination method depend 

on the light intensity [5]; the values of R determined 

s 

by this method fall in the range indicated in Table I. 

In contrast to the illumination method, method (a) involves 
only electrical eKcitation and applies to ordinary diodes 
as well as to solar cells. 

The high-frequency measurements were done on Wayne-Kerr 
B601 and B801 bridges; low-frequency measurements were done 
on a Wayne-Kerr B224 bridge. 

3.4 Summary 

This work demonstrates a simple and rapid method 

(method (a) ) for determining the series resistance of diodes 

and solar cells based on the measurement of small-signal 

admittance. The OCVD method (method (c)) provides a rapid 

but coarse estimate of R . Method (b) is useful for devices 

s 

having a short substrate lifetime. Method (b) enables also 
measurement of R^ at any bias voltage V, which allows the 
determination of any dependence that R might have on V. 

5 
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Series 

resistance R 

s 

in ohms 

Device 

Method (a) 

Method (c) 

Photovoltaic 
output [5] 

+ 

p -n diode 




lOficm 

20 

19-28 

15-23 

p"*'-n diode 




0 . If2cin 

4.35 

3-10 

3. 5-4. 5 

2 + 

2x2 cm n -p solar cell 


• 


2Sicra 

0.33 

-0.6 

0,25-0.34 


Table I 


Series resistance from methods (a) and (c) together with results 
using light excitation for comparison. 
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Figure 1 (a) ac equivalent circuit; (b) high-frequency equivalent circuit for wCj) >> Gj^; 

(c) parallel equivalent circuit as measured by bridge, where Gp, Cp are functions 
of Go, Co, Rg. w; (d) dc equivalent circuit. 
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CHAPTER 4 FORWARD-BIAS CAPACITANCE AND CURRENT MEASUREMENTS FOR 

DETERMINING LIFETIMES AND BAND NARROWING IN P-N JUNCTION 
SOLAR CELLS 

4.1 Introduction 

One of the material parameters that most influences the performance 
of a p-n junction solar cell or diode is the minority-carrier diffusion 
length L of the base (or substrate) region. The diffusion length limits 
the collection efficiency and power conversion efficiency achievable in 
a solar cell, and it controls the current-voltage characteristic in a 
long-base (nonlUuminated) diode. 

Many methods exist for the measurement of L, Of these, the ones of 
most relevance to p-n junction device studies yield the value of L in a 
solar cell or diode after the junction is formed so that the influence 
on L of the process used to fabricate the junction is included. Existing 
measurement methods applicable to this problem [l]-[5] involve various 
tradeoffs between accuracy and convenience. 

The purpose here is to describe a new method [6] applicable to solar 
cells for measuring L and thus the minority-carrier lifetime T, which is 
accurate and requires only straightforward capacitance measurements made 
at the device terminals. The accuracy results because the measurements 
required employ very accurate capacitance bridges. The accuracy is demonstrated 
by showing that the values of L yielded by this capacitance method agree 
well with those obtained by use of the X-ray method described by Rosen- 
zweig [1], In contrast, the results determined by the capacitance method 
disagree sharply with those determined from the open-cir cult-voltage-decay 
method [4], for reasons that are discussed. 
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In its present form, the irethod measures L for a non-11 laminated 
cell. We discuss the conditions under which this value of L is appropriate 
to use in describing an illuminated solar cell, and indicate experiments 
that confirm these conditions. 

The method can also determine material properties of the quasi-neutral 
emitter region. We discuss the determination of the energy-band narrowing 
and the phenomenological lifetime in the emitter both by this method and 
by another method involving only dc measurements of emitter current. 

4.2 Theoretical Grounds for the Method 

Consider a p -n junction diode or solar cell with no illumination 
applied; analogous results hold for n -p devices. The minority-carrier 


diffusion length, L = /d T , in the n-type base influences the small- 

It y It 

signal capacitance and conductance Ggjjjg of the quasi-neutral base 

region according to [7] 


Aqn?L 

‘^QNB "" M [expC^) - 1] 


Cl) 


AqnTD „ 

= A. -R [exp(-3^) - 11 


QNB kT 


( 2 ) 


for an applied sinusoidal signal whose period greatly exceeds the base 

lifetime T . The above results hold for a one^dimensioanl model of a 
P 

device having a uniformly doped base and operating under low-injection 
conditions. The undefined symbols appearing in (1) and (2) have their 
usual meanings. 

If either or could be determined from measurements, (1) 

and (2) would then yield L^, which is the desired quantity. But the 
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capacitance C and conductance G measured at the terminals contain com- 
ponents in addition to C„_ and G„„_ [8]- [11] 


C = 


G = 


Cl + Cg + Cgg + 


G_ + G__ + G„„ + 

S SC QNE QNB 


(3) 

(4) 


under low-injection conditions. In (3) and (4) the subscript S refers 
to the surface region, SC to the junction space-charge region, and QNE 
and QNB to the quasi-neutral emitter and base regions, respectively. 

The capacitance C^ is the capacitance due to the ionized impurities in 
the junction space-charge region and Cg^ is that due to the mobile holes 
and electrons in that region [8]. From the viewpoint expressed in (l)-(4), 
the problem of determining becomes the problem of extracting from 

the measured C versus V characteristic or of extracting from the 

measured G versus V characteristic. 

As experiments [12], [13] including those described in Section III 
show, G can be comparable to G_™ or larger than G in many p-n 

QNE QJNiJ QNE 

junction devices of interest. Thus, except in special cases, one cannot 

determine G_„_, and hence, L , from the conductance data. 

QNB ’ p 

The key to the extraction of C , and hence to the determination 

QNB 

of L , from the C-V data lies in making capadtance measurements at 
relatively high as well as at low frequencies. At low frequencies, the 
mobile carriers in all regions of the device respond to the signal, and 
(3) describes the measured low-frequency capacitance. But, at a higher 
frequency only mobile carriers associated with Cj and Cg will follow 
the signal [7] and will contribute to the measured capacitance if the 
period of the signal is long compared with both the dielectric relaxation 
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time of tfie semiconductor and the transit time across the junction space- 

charge region. Thus by subtracting this higher-frequency capacitance 

from the low-frequency capacitance, one obtains a capacitance associated 

with the quasi-neutral regions. This capacitance is if the 

phenomenological lifetime Tg of the emitter [12]- [14] is much shorter 

than Tp, and if the signal frequency f is chosen so that Tg«l/f«Tp. 

If, on the other hand, T„-T , or if l/f«T_,T , this subtraction yields 

Up E p 

the capacitance, + ^QNB’ reduces to for the commonly 

occurring case that Thus as will be discussed and illustrated 

QNE QNB 

in the following sections, the method yields C „ and, hence, from (1), the 

QNB 

value of L , 

P 

The method applies directly to large-area solar cells. In soma cases, 
however, advantages may result from making the measurements on sraall-area 
test devices. For example, this may be desirable if high-frequency ad- 
mittance bridges having a suitable capacitance range are unavailable in 

a particular laboratory. For small-area devices, with radius a~L , the 

P 

two-dimensional fringing effects illustrated in Fig. 1 invalidate the 
one-dimensional theory underlying (1) and (2) and corrections must be 
inade. Correction curves were calculated [15] by solving the continuity 
equation for holes in tx^o space dimensions. One set of correction curves 
is shown in Fig. 2. 

The theory underlying the method, contained in (1) and (2), applies 
to non-illuminated solar cells. Solar Illumination can change the carrier 
lifetime because the solar spectrum can change the transitional probability 
rates of the intermediate levels in the forbidden band [16]. In general, 
the method yields values of lifetime appropriate to an illuminated cell if 
the shifting approximation holds, that is, if the current under illumlna- 
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tion is the current in the dark shifted by the short-circuit current [17]. 
Very careful measurements made on various single-crystal silicon p-n 
junction solar cells demonstrated the applicability of the shifting 
approximation for AMO and lower levels of illumination, thus indicating 
the validity of the method for these cells at these illumination levels. 

In assessing the applicability of the shifting approximation, corrections 
may be needed to account for the effects of series resistance, 

A quick check for the light- dependence of the lifetime can be made 
by measuring the forward-biased capacitance in the dark and again xrtth 
the light on. 

4.3 Illustrative Example 

To explore the utility of the method, we fabricated s set of small 

circular p -n junction solar cells having various base doping concentra- 

3 17 3 

tlons between 1.4x10 cm and 6.5x10 cm . Each diode x?as circular 

-3 2 

with an area of 4.5x10 cm . Three separate runs were made; runs 1 and 

3 at 950°C and run 2 at 1000°C. Runs 1 and 2 employed phosphorus gettering 

from the back face to getter metallic impurities from the substrate. 

19 —3 19 —3 

Boron surface concentration ranged from 1.5x10 cm to 4x10 cm 

and the junction depth ranged from 0.3uni to 0.9pm. 

Figure 3 shows representative data: the measured f»n(G) and 2.n(C) 

dependencies on qV/kT for device from run 1 having a base doping concentra- 
15 -3 

tion Npjj = 1.25x10 cm . The data x^ere taken at 299 K and 376 K. The 

low-frequency capacitance and conductance, and x^ere measured at 

3 kHz using a Wayne-Kerr B224 bridge. High-frequency capacitance was 

Hr 

measured at 10 MHz using a Wayne-Kerr B801B bridge. The level of the ac 
signal was kept below 5 mV. 
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To calculate the diffusion length L of the base region, we assume 

P 

that the minority carrier charge storage in the base greatly exceeds that in 
the emitter, and, hence, that ^or the particular diode under 

study, these assumptions are likely to hold for three reasons i a) the 
low surface concentration Kg tends to limit the size of the energy-gap 
shrinkage effects [12], [18], [19] in the emitter; b) the low doping con- 
centration Njjjj of the base Implies a large minority carrier charge storage 
in the base; c) the thinness of the emitter region tends to limit the 
amount of minority-carrier charge it can store, 

Cqjj = is obtained by subtracting from which yields for 

iln versus qV/kT a straight line having unity reciprocal slope. Then, 
from (1), the minority carrier diffusion length of the base obtained 

from the capacitance method is L „ = 80iim. To assess the accuracy of this 

pc 

result, we also determined the diffusion length by the X-ray method 

described by Rosenzweig [1], which is accurate within about ±5 percent [1], 

and found Lpj^ = 84iam. The good agreement between and Lp^^ tends to 

confirm the utility and accuracy of the capacitance method. 

For some devices, the emitter recombination current wi3.1 be much less 

than the base recombination current, so that and L can be 

QNE QNB p 

deduced from (2) using only the low-frequency conductance data. This is 

true for the device under consideration in Fig. 3. If one assumes 

one obtains L „ = 80]im, which agrees well with the values 
QNE QNB pG 

obtained by the capacitance and X-ray methods. was obtained using a 

construction technique described elsewhere [13]. 

But the conductance method falls whenever the emitter recombination 
current contributes significantly to the current-voltage characteristic, 
as can be expected to occur in diodes having higher base dopJ ng c(7n- 
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centratlons and consequently lower minority-carrier charge storage in 
the base. This point is illustrated in Fig. 4. For diodes having base 
doping concentrations above 10 cm , differs from and by 
a factor of two or more. The capcitance and the X-ray methods give 
results in excellent agreement for all base doping concentrations, how- 
ever, which Illustrates the accuracy and applicability of the capacitance 
method. Excellent agreement between Lpg and ^PX also justifies the use 
of dark Cgjjg for determination of Lp in solar cells, since X-ray radiation 
can excite traps in the bandgap in a similar way as light excitation - 

For our devices the ratio a/Lp is larger than about 4 and the two- 
dimensional correction factor is small, as can be seen from Fig. 2. 

Appropriate corrections were made, 

A. Comparison with the Open-Circuit-Voltage-Decay Method 

Determination of the base lifetime T by the open-circuit-voltage- 

P 

decay (OCVD) method [4] is a commonly used approach. The OCVD method 
is attractive because it involves the measurement of only a single 
transient, which can be rapidly done. 

To examine the applicability of OCVD to silicon devices, we consider first 
its theoretical basis. In the theoretical development of the OCVD response, 
one finds the open-circuit response of a diode following the sudden removal 
of forward bias by solving the continuity equation for the minority carriers 
in the base. Under low-injection conditions, the open-circuit voltage 
decays linearly with time, to a close approximation, and measurement of 
the slope of this straight-line decay suffices to determine This result 

applies if the minority carriers in the base dominate in determining the 
response, as was the case for the germanium diodes considered in the 
original OCVD study [4]. But if i.he carriers associated with the junction 
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space-charge region or with the quasi-neutral emitter contribute 
appreciably to the response, then the OCVD method is not applicable. 

The junction space-charge region contributes to the measured 
capacitance and conductance characteristics of Fig. 3, as evidenced 
by the nonunity slopes of these characteristics, particularly of the 
capacitance characteristic. At 299 K, the slopes come nearest to 
unity in the voltage range, 0.45 V < V < 0.52 V. For this range, the 
minority hole concentration in the base is in low injection. Figure 5 
displays the corresponding OCVD response, taken at 299 K and beginning 
at V = 0.5 V. Note that this response fails to exhibit a simple linear 
decay. Thus as one could anticipate from the nonideality of conductance 
and capacitance characteristics, the OCVD method for determining is 
inapplicable for this diode at 299 K, A course estimate of can be 
Inferred from the OCVD response, however, by fitting a straight line 
to the response and measuring its slope. This yields x^ (OCVD) ^ 8|is, 
which differs by nearly a factor of two from the value, x^ - 4.4ps, 
corresponding to = 80ym determined by the capacitance method and 
corroborated by the X-ray method and by the calculation [7], x^ = . 

For comparison, the straight-line decay corresponding to Xp = 4.4ys is 
plotted on Fig. 5. 

According to theory, the quasi-neutral capacitance and conductance 
increase more rapidly with temperature than do the components from the 
junction space-charge region. Thus one might expect the capacitance 
and conductance characteristics to have slopes that more closely approach 
unity as the temperature is raised. This is confirmed by the data of 
Fig. 3 measured at 376K, where the closest approach to unity slope occurs 
in the voltage range, 0.25V < V < 0.32V. In this range, which corrospondH 
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to low-injectioa levels of the minority holes in the base* the conductance 

almost shows an ideal slope of unity* and the capacitance shows a 

reciprocal slope of 1.07. As Fig. 5 indicates* however, even these 

slight departures from the ideal characteristics are enough to produce 

an oeVD response that again fails to ejdiiblt a simple linear decay. 

The fitting of a straight line to the OCVD response at 376 K results in 

the estimate X (OCVD) - 9{Js, which differs considerably from the value 

X^ a 7ys, corresponding to “ SOpm at 376K determined by the 

capacitance method and corroborated by calculation [7], x^ = ^^qN^^QN' 

Similar deviations from the ideal capacitance and conductance 

characteristics and the ideal (linear) OCVD response were seen for all 

of the silicon diodes of Fig, 4. This suggests that the OCVD method 

for determining x is questionable for silicon devices, although one 
P 

might expect its accuracy to increase as temperature is raised. The 
capacitance method avoids the inaccuracies resulting from the nonideality 
of the characteristics by combining low-frequency and high-frequency data. 

As was discussed previously [13]* the junction current recovery 
method [20] is generally even a poorer method than OCVD for determining 
carrier lifetime and diffusion length of the base region. 

4.4 Discussion of the Capacitance Method 

The capacitance method for determining and x^ involves measuring 
the dependence of the capacitance on forward bias at two frequencies: 
a frequency fj^, low enough that and a frequency fj^ high enough 

that l/^H^Xp, Once the low-frequency and higher- frequency capacitance 
are measured, they are subtracted to give the quasi-neutral capacitance 
Cqjj. The resulting dependence of JlnCC^^) on qV/kT must be a straight 
line of unity slope, which provides a convenient self-consistency check 


piVGE IS 



-43- 


on the validity of the results, For diodes for which which 

often occurs, the quasi-neutral base capacitance C„„_ is thus determined. 

Combining C„„„ with (1) then determines L and t . 

QNB p p 

The inequality, for the cells of Fig, 4 because 

19 —3 

the surface doping concentration is low (Hg - 10 cm ) and the emitter 
junction is thin (x^ ~ 1 ym). These conditions make the total number of 
dopant atoms in the emitter small. That Cqj]g«CQjjg for these devices 
is substantiated by the demonstration that L „ = L But if a given 

pij pX 

device has a large number of dopant atoms in the emitter, then can 

be large because of band-narrowing effects [12], [18], [19], If, in 
addition, the device has a large doping concentration in the base, which 
from (1) implies a small value of C_„_, then = C_„_. For such a 

device, the diffusion length of the base can be determined by the X-ray 
method, and the method described in this paper can be adapted to enable 
the determination of the minority-carrier charge Qg stored in the emitter 
and the phenomenological lifetime of the emitter. From the values of 
Qg and Tg thus determined one can deduce information about the band- 
narrowing and recombination mechanisms in the highly doped emitter [12], 
The details of such studies for bipolar transistors were published else- 
where [21], [22]. Similar work on solar cells is in progress. 

There are several practical limitations that define the scope of 
applicability of the capacitance method and they are discussed in detail 
in Ref. [6], These limitations are not severe and the capacitance method 
applies to a wide range of silicon devices, as the results given in 
Section III demonstrate. The capacitance method has the advantage of 
being accurate and of requiring only straightforward electrical measure- 
ments at the terminal using Inexpensive equipment. It also enables the 
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easy determination of the dependence of diffusion length and lifetime 
on temperature. 

The accuracy of low-frequency capacitance bridges, such as Wayne-Kerr 

B224, is about ±0.1 percent and the accuracy of available high-frequency 

bridges is about ±2 percent. If the temperature during measurement is 

2 

controlled within ±0.1 K, which will introduce an error in n^ less than 
2 percent, the overall accuracy of the method is estimated to be about 
±5 percent. 

4.5 Bandgap Warrowing in Emitter from Temperature Dependence of Emitter 
Current 

Slotboom and deGraaff [23] and Martinelli [24] recently proposed 

methods to determine bandgap narrowing AEg in the base region of bipolar 

transistors. Because the base doping of transistors is limited to about 
19 -3 

10 cm , these methods do not allow investigation of an entire range 

21 -3 

of high doping densities up to about 10 cm . Although these and similar 
methods provide useful information about bulk material remote from a sur- 
face, they cannot directly yield AE^ of the thin emitter region of any 
particular device under study because AEg may depend in part on the prox- 
imity of the semiconductor surface near which lattice strain and impurity 
and defect clustering can occur. In addition, AEg may also depend on the 
particular diffusion or ion-implantation processing used to form the 
emitter region. 

To determine AEg of the emitter, we propose an alternate method 
which employs the temperature dependence of the emitter current. This 
method applies over the entire range of emitter doping concentrations 
present in p-n junction devices. If the emitter doping profile is flat, 
the interpretation of the results of the method is particularly straight- 


forward . 
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The simplest version of the method results if the emitter is 
transparent, i.e., if Injected holes (for an n -p device) recombine 
mainly at the surface rather than in the bulk of the emitter. If 
the emitter region is thin, the transparency of the emitter can be 
demonstrated experimentally even for very high emitter doping con- 
centrations. A theoretical and experimental discussion of this issue 
will appear elsewhere. 

If the emitter is transparent and if the surface recombination 
velocity Sp is large, then [25] 


r W. 




E 


t~l 


N(x) 


dx 


0 D (x)n. (N) 
p le 


(5) 


2 3 

where exp(-EgQ(N)/kT) in which C is a constant nearly independent 

of the doping concentration [23]. EqqCN) is the bandgap extrapolated to 

0 K, which equals 1.206 eV for low doping concentrations [23,26]. For 

a degenerately-doped emitter region, in which the electron concentration 

N is large, the mobility is assumed to be independent of temperature. 

2 

The quasineutral emitter thickness is W^. In (5), n^^ is the effective 
intrinsic density. 

For a flat doping profile, (5) becomes 


- =i^'‘exp(-l!go(N)/kT) 


( 6 ) 


where is independent of temperature. 

Thus a plot of fi-nCJ^^/T*^) versus 1/T yields a straight line from which 

EO 

can be calculated. The emitter saturation current is deter- 
mined from the measured total saturation current by teHinlf|ut‘8 rtcHi* r Lbud 
elsewhere [12,13] . 
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The method was used to measure E_ on a device with almost flat 

b 

20 -3 

implanted As profile with N - 1.5 x 10 cm and - 0.37pm. The 

plot of Xrn(J_„/T^) versus 1/T for this device appears in Fig. 6. From 
£0 

the slope = 1.046 eV which gives = 1.206 eV - 1.046 eV = 0.160 eV. 

17 

Additional points for the entire range of emitter dopings from 10 
-3 21 -3 

cm to 10 cm are being measured and adjusted to . include the effects 
of Fermi-Dirac statistics. The results will be published elsewhere. 
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Fig. 1 Cross-section of a diffused p -p junction diode showing a 
two-dimensional spreading effect. Boundary conditions for 
hole concentrations at the junction edge P(x.) and at the 
back contact P(W) are indicated. W is the base width and 

S is a surface recombination velocity at the back contact. 
P 






ni = 1 .01 


!<: — m=i . 11 


tn = 1. 0 




- '"hf 


V (Volts) 


m=:1.33“^ 


'Tn=1.0 


V (Volts) 


Measured conductance and capacitance versus forward bias V 
for diode with = 1.25xl0l5 cm"3. a T=376 K; b T=299K. 
The insert shows the capacitance characteristic for diode 
with Njjjj=8x 10^6 cm“3, demonstrating the increased importance 
of Cgji and emphasizing the value of the subtraction for 
determining Cgjjg. In the figure, — stands for the com- 
ponents from the quasi-neutral regions (m=l) and 
stands for G + G„. 









0.6 



Time (pser) 

Fig. 5 OCVD response for device of Fig. 3 at 299 K and 376 K. Dashed lines correspond to 
the lifetime T determined by the capacitance method. 
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CHAPTER 5 

AN ANALYSIS OF THE DIFFUSION CURRENT AND CAPACITANCE IN A TWO 
DIMENSIONAL p-n JUNCTION DIODE FOR DETERMINING 
THE DIFFUSION LENGTH 


5,1 Introduction 

In, a previous related publication [1] by the 
authors a method was proposed for determining the diffusion length in 
the base of a p-n junction diode from measurements of the diffusion 
current-voltage and diffusion capacitance-voltage characteristics. As 
a central approximation in the development of this method, the current 
and capacitance characteristics were assumed to be those resulting from 
the conventional one-dimensional model of the diode, in which excess 
hole and electron distributions are adequately described by a single 
space variable x, the distance from the junction surface. This approxi- 
mation describes the characteristics with sufficient accuracy if the 
area is large, or more precisely, for a circular device, if the radius 
a of the emitter region greatly exceeds the minority-carrier diffusion 
length Lp of the base region. But the bridge used for the high-frequency 
capacitance measurements required in this method imposes an 
upper limit on the capacitance that can be measured. Thus diode 
structures must sometimes be used for which the approximation is 
liivulLd and for which the one-dimensional model of the diode intro- 
duces Intolerable error in the detenninatlon of L . Figure 5.1 

P 
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lllustrates a p-n junction diode that* if were large enough* would 
be poorly duscrtheU by n one-illinanaimml model. 

The purpose of this chapter Is to provide .ui auftlyalK of tho dif- 
fusion current and capacitance characteristics for a two-dimensional 

model of the p-n junction diode that allows the determination of L for 

P 

any value of a/L . 

P 

Grimbergen [ 2 ] has treated a problem related to the one of interest. 

Grimbergen studied a p-n-o structure. He assumed the surface recomblna- 

4* 

tion velocity at the high-low n-n junction to be zero and the base 

width W to be much smaller than L^. Grimbergen reduced the two- 

dimensional problem to two coupled one-dimensional problems. As will 

be shown, the result he obtained is close to one of the cases of our 

two-dimensional calculation, that corresponding to = 0, and 
W 

r — - 0.1 in Fig. 5.2. But, for many diodes, the assumption S = 0 at 

J-i P 

P 

the back, face is inappropriate and the ratio W/L^ is not necessarily 
small. Thus, a general solution of two-dimensional problems will be 
useful. 

The chapter begins by defining the boundary-value problem. It 
resembles boundary-value problems encountered in the conduction of heat 
but differs in that a volume recombination term appears in the dif- 
ferential equation. The solution of this boundary-value problem is then 
presented, which yields the current 1(a) and the capacitance C(a) 
normalized by their one-dimensional counterparts plotted against the 
ratio, a = in conjunction with current and capacitance 

measurements, the plot enables the determination of L . This is illus- 

P 

trated for two devices with different areas. The diffusion length L 

P 

determined by this method is shown to be consistent with that determined 
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by the X-ray method [ 3j which illustrates the good accuracy to be 
expected from use of the two-dimensional analysis. 


5.2 Theoretical Analysis 


The hole diffusion problem in the n-type substrate of a shallow 
diffused p-n junction (see Figl 5,1) is similar to the heat conducting 
problem of a metal disk [4,5 j, except that a voltime recombination terra 
appears in the differential equation, describing the kinetics of the 
holes. The differential equation originates from the continuity equa- 
tion, which for the steady-state condition is 


with 


V-J 

P 


-qP'/T 


P 


J " -qD VP 
P P 


Combining (5.1) and (5.2) yields 


(5.1) 


(5.2) 


7 

V“P 



(5.3) 


where P’ = P - P is the excess hole concentration and L = >(D t is 
no p p p 

the diffusion length of holes in the base region. In cylindrical co- 
ordinate (p,(J),z), equation (5.3) becomes 





1 ^ 
P 3p 






(5.4) 


Separation of variables is accomplished by substituting P*(p,(j), 2 ) = 
R(p)4>(d))ZCz), which leads to 
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1 0 + M 

R 3^2 Rp 3p 


i -L 
2 
P 


2 2 

3 <!> 1 3 2 

2 Z 2 

3<f 3z 



0 


(5.5) 


Normalization with respect to the junction radius a can be done by intro- 
ducing p - ra. Thus, equation (5.5) becomes 


2 2 2 2 

i ill + X M 4. 1 XXI + a_ ±_z 


Now, let 


and 


ai ifz a^ , 2 

2 3 z 2 " L 2 ■ ^ 

P 


1 2 
0 “2 " 


= a 


Substituting (5.7) and (5.8) into (5.6), we obtain 


X_| + i ^ + (fc2 _ V ^ 0 


(5.6) 


(5.7) 


(3.8) 


(5.9) 


(5.10) 


For a problem possessing azimuthal symmetry, u = 0, i.e., <& = constant. 
Thus, (5.6) reduces to the two ordinary differential equations (5.7) 
ard (5.10). Solutions of (5.7) and (5.10) are 

/JT~2 /. 2 ^ 2 

/k+a Y k + a 

- z — z 

Z=Ae ^ +Be ^ (5.11) 


R = CJQ(kr) 


(5,12) 


rher afore, we have 


P'(r.z) 


/. 2 ^ 2 /, 2 ” 2 

vk+a vk+a 

to a> _ — — 2 z 

J R*Z - e "" + D e ) J.{kr) 

k»0 k*0 


(5.13) 


For a diode with 
the back contact 


base width W and a surface recombination velocity S at 

P 

(x = W) , the boundary condition at z = W requires that 


D 

P 


3P* (r,z) 
dz 


z=W 


-Sp P'(r,z) 


z=W 


(5.14) 


The negative sign in (5.14) is defined by the current direction. Sub- 
stituting (5.13) into (5.14), for any number of k, we have 


D [(-A) 
P 


42. 2 - 




k -ta* 


+B 


/k^-Kt^ 


k +a" 


W 


/TU 

■ / k +a 


*1 

J = 


-S [a e 
P 


W 


/? 


k +a 




+B e 
(5.15) 


From (5,15) we obtain 


B = 


D / 2 , 2 

- P , / k + a 
1 - c * 

S a 

P 


/k^ 


2 W 


D > 72 ” 2 

, , p k + a 

^ + s — 


(-A) 


(5.16) 


Substituting (5.16) into (5.13), we have 
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P'Cr.z) = I A[e 
k=0 


.'52, 

® ] Jrt(kr) 


z S a 
3 g— 


i+j — ^ 
p 


.^^ 2„52 


= /A(k)[e 
0 


1+ -JP , +a 
S a 

P 


3Jq ( kr)dk 


(5.17) 


For a shallow diffused junction diode » with zero surface recombination 
velocity at the front side of Si-SiO^ interface, the boundary conditions 


require that 


S a 

P'(r,0)=P'(x.) = / A(k)Cl B 


1- ^.52 


2W 


° /TT2 

1+ 

for 0 < r < 1 


]jQ<kr)dk (5.18) 


and 


°o = ^ 




2W 


D /T “ 2 

1+ ^ 

S a 

P 

for 1 < r < ” 


]j^(kr)dk 


(5.19) 


Such pairs of integral equations are known as "Dual Integral Equations" 

[6], We now compare (5.18) and (5.19) with the following integral pairs: 
00 

/ G(k)f(k)j.(kr) - Ar*^ - A 0 < r < 1 
0 ^ 


(5.20) 
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and 


/ f(k)J^(kr) - 0 
0 ^ 


1 < r < 


Thus 


f(k) = A(k) 




lA.tlZZZ / 


S a 

1 + E 


, 2 ^ 2 
k + a 


2W 


1+ Vk^ + o.^ 

% ^ 


and 


G(k) ^ 


/^T7^ 


1-- E. 


2W 


■ 


/. 2 . 2 
/ k + a 




1 +- 


2W 


If we choose 


"p 


f(k> = k^'2 I ,(k) 


m- 


Q m 2mHJs' 


determined from the following simultaneous equations J 

CO I — “ 

Jo ^m-o\ "ill for n = 0 


a + L a =» 0 
n m*n n 


for n > 0 


where 


Numerical methods are used iu computing L 

m*n 


dk 


(5.21) 


(5.22) 


(5.23) 


(5.24) 


(5.25) 


(5.26) 
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Once is determined, the total hole diffusion current across the 
junction is 


I = qD / 2irpdo “ qD 2irrdr 

P 0 P 0 3^ 

After integrating (5.27) and dividing it by the one-dimensional current 
Iq, which can be expressed as 


qD „ 

Iq " * Tra P"(x ) 

P ^ 


(■■■ P P -)cosh(^) + sinh(^) 
2 E E_ 

(- P P -)sinh(~) + cosh(“) 


(5,28) 


P 


the ratio of I /I^ appears as plotted as a function of a in Fig. 5.2. In 

CL U 

this figure, W/L and S L /D are parameters. 

P P P P 

For a long-base diode, i.e., one for which W >> L , the constant B 

P 

in (5,16) approaches aero. For this case, the detailed solution of the 
two-dimensional boundary value problem is given in Appendix II of Ref. [7], 
result is [8,9] 


f (a) 


I 

a 




(Aq + O.Aj^ + O.A2 + 


. - ) 



(5.29) 


Because the series Z A converges fast, neglecting terms higher than 

m=0 ® 

yields accuracy of i3%. The same procedures can be 

applied to calculate f(a) for different values of W/L and S L /D . 
t'i' P P P P 


The 



Figure 3.2.1 


Normalized current density factor f(a) = Jq/Jo’ ^ 
function of the ratio of junction-area radius and dif 
fusion le igth for Sp = « at back contact. 


s 


f (ct) calculated for -z — ■ 0,1 

£j 

D P 

r?" ■ " ®p ■ *» 

n n * 


S L ^ 
P P 


C: -r-f- = 0 (i.e. S = «>) 

b L D 

p p 


Grimbergen's approximation 

(with S- - 0, W-x./W - 1, 

W <Lp)^ 


Figure 5.2.2 Illustrating that the detailed solution presented here 

(curve A) agrees well with Grimbergen's approximation tor 
diodes for which W/L„ << 1, S_ = 0, and’ (W-x.)/W = 1. 
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5,3 Experimental Results 

Using the exf urlmental set-tip of Fig, 5. 3, we measured the current- 

Voltage characteristic of p -n junction diodes with different areas, described 

in Table 5.2, For the measurement, the devices were heated inside a 

Statham furnace to a high enough temperature that the base recombination 

current dominates that from the junction space charge region, as 

evidenced by a unity slope on the (In I) vs. (qV/kT) characteristics. 

The existence of a unity slope indicates that the current I is the sum 

of the recombination current from the quasi-neutral base and emitter 

regions. From the experiments [10] and reference [ 1 ], we know 

that the base recombination current dominates for our devices having a 

16 “3 

base doping concentration less than 10 cm . 

Table 5.1 shows the current measurement on devices with seven dif- 
ferent areas. Measurements are taken at a temperature of 375.5 ± 0.1°K 
and at a fixed forward bias of 0.3 volts. The base diffusion length of 

the devices obtained from X-ray measurement is L =80 um, and the sub- 

px 

15 -3 

strate doping concentration is 1.25 x 10 cm 

Experimental results are compared with the theoretical predictions 
in Fig. 5.2. Good agreement is obtained. 


5.4 Determination of Diffusion Length 

As Chapter 3 has discussed, the diffusion length can be determined 

by (a) current measurement or (b) capacitance measurement. The methods 

of Chapter 3 hold directly provided ~ >> 1, i.e., they hold for the 

P 

one-dimensional case only. If the device is small enough that — >> 1 

is not true, then methods (a) and (b) must be altered as follows to 

include the dependence of P and J m two space dimensions. 

n p 


PAGE IS 

quality 
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Figure 5.3 Block diagram for current measureraeat 


Table 5.2 


Diamecers 
in Mills 


I 

a 


a 



tra 

J 


a 



f(«) 


Diffusion current measured for seven devices having different values of ct = a /Lp^,. 
For these devices the one-dimensionul current density is 

n 2 HY 

i IfT -T 

J, = qD — — (e - 1) = 5.91 K 10 amp. 

The measured correction factorj shows excellent agreement with the correction 

fac tor, f (ct) , calculated from the solution to the two-dimensional boundary-value problem. 
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30 

20 
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7 
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0.30x10“^ 
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1.58 

1.11 

7.02x10“^ 

7.89x10"^ 

9,27x10^^ 

1.21xl0~^ 

1.19 

1.34 

1,57 

2.04 

1.22 

1.36 

1.75 
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Method (a): For devices for which the base recombination current 

greatly exceeds the emitter recombination current (for our devices if 
16 "3 

Ng i 10 cm ), the diffusion length can be determined by the following 
scheme. Suppose currents are measured (see Section 3) on devices 1 and 
2 of radius a^^ and a^, then 


2 

I ^ A.J = ira.J ^ A.^ J_ = A.J. = A.J„£(a.) 
j J “j J ° ^ ° J 0 J 


a . 


(j “1 or 2) 


(5.30) 


where a. - and 
2 ^ 

p 


n. 




T i 1 kT C V 

'-“8a) 

DD p P 


(5.31) 


There are two equations and three unknowns. To find the value of an 

iterative procedure is used as follows. Try an initial value of L and 

^2 ^ 
calculate U 2 = find the fCa^) from the curve of Fig. 5.2. Plot 

f(a, ) = f(a„) Notice whether f(a.) is on the curve or not. If 

not, choose another L , check f (a, ) again. Continue this procedure 

P 1 

until f(a^) and ^( 02 ) tioth lie on the curve. The value of Lp producing 

this condition is the correct value of Lp. Note that good accuracy 

results provided and U 2 differ considerably: For example, ~ 

and a, = 1 yields accurate determination of L . 

1 ^ P 

Method (b): If the emitter recombination current is comparable to 

16 ***3 

that of the base (for N„ 10 cm ), the above scheme is inappliacable. 

D 

However, the capacitance method [1] still 


can be used to determine L . This method determines the base diffusion 

P 

capacitance [ ll] which, for devices 1 and 2, is 
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«D 


dv 


dl 


D. 

1 


dV 


— -3- j 
T, kT 

J J 


(J - 1 or 2) 


(5.32) 


But devices fabricated on the same wafer. Thus the deter- 

mination of L from C_. and C_„ is essentially the same as the deter- 
p ul UZ 

mination from 1^^^ and described in (a) above. The theoretical curve 
of Fig. 5.2 therefore can also be used for the capacitance method of 
determining L^. 

In contrast to the methods of Section 5. A, which require two or more 
devices in determining the diffusion length, the diffusion length can be 
determined from measurements made on a single device by use of the 
following iterative procedure. First, is estimated and the cor- 
responding one-dimensional capacitance (or current) is calculated from 
(3.11) and multiplied by f(a) from Figures 5. 2. 1-5. 2. 3. The result is 
compared with the measured (two-dimensional) capacitance (or current). 

The procedure is continued until self consistency results. 


5,5 Discussion 

The diffusion length as a function of temperature T can also be 
determined by use of Fig. 5.2, In contrast, the X-ray method does not 
readily yield as a function of T because the small short circuit cur- 
rent (due to the small generation rate of X-ray) is obscured by the 
thermal noise signals existing in the device and contact wires. Thus 
the X-ray measurement of is difficult at temperatures higher than 


room temperature. 
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Determination of the diffusion length by the current measurement 
is simple and accurate. 

The validity of the theoretical calculation presented here can be 
checked by letting a 0, and considering the special case in which the 
volume recombination vanishes. This produces the electric disk problem 
[ 53 , illustrated in Fig. 5.1. For ci = 0, = 1, =* . * 0. 

Now let the earth ground plate be placed near the electric disk. For 
example, let a/h = 20, where a is the radius of the disk and h is the 
distance between disk and earth plate. Then, the average electric field 
calculated using section 5.2 is 

E =‘~x 1.08 (5.33) 

ave h 

where V is the potential of the disk. But, if a/h >> 1, we know from 
electrostatics that E = V/h. Thus calculations using Section 5,2 are 
shown to agree with the results of the analogous electric disk problem. 

In our theoretical calculation, we have approximated the shallow 
discussed junction by an infinitesimally thin disk. The current density 
becomes infinite at the edge of the junction, which is contrary to the 
condition occurring in actual diodes. The junction depth shown in 
Fig. 5.4 is small but not equal to zero. As is shown in Appendix VI of Ref. [7], 
the current density at the junction edge can be expressed by 

X. 

■^EDGE “ 

where and are the modified Bessel functions of the third kind with 
orders zero and one, respectively. 
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INCREASING Vqq BY THE HLE STRUCTURE 


PART III! 



CHAPTER 6 EMITTER CURRENT SUPPRESSION IN A HIGH-LOW-JUNCTION 
EMITTER SOLAR CELL USING AN OXIDE-CHARGE-INDUCED 
ELECTRON ACCUMULATION LAYER 
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INTRODUCTION: As Brandhorst first noted [1], the power oonvevslini ofi'icloiic'y 

ri seen In silicon p-n junction solar cells is considerably less than the 

maximum theoretical value of n mainly because the open-circuit voltaRC 

is smaller than simple p-n junction theory predicts. Experiments [2] on 

n -p silicon cells have shown that this discrepancy in results from the 

dominance, in the non- illuminated (dark) cell, of the emitter recombination 

current Jjj over the base recombination current J^. In cells having base 

1 / "“3 

doping concentrations of the order of 10 cm , for which the largest 

values of V_„ are seen, J_ exceeds by about an order of magnitude [2], 

rather than being several orders of magnitude less than as Is predicted 

by simple p-n junction theory. The excess Jg has been attributed to the 

mechanisms [3] of energy band-gap narrowing and lifetime degradation that 

accompany heavy doping concentrations in the n emitter. 

To suppress and thus raise the achievable value of a new 

structure, the HLE junction solar cell, containing a high-low (H-L) junction 

in the emitter, has been proposed and its performance has been calculated 

on theoretical grounds [4,5]. The purpose of the present paper is to 

demonstrate, by experiment, that an H-L junction in the emitter created by 

an oxide- charge-induced electron accumulation layer can completely suppress 

Jg. We call this device the oxide- charge- induced, high- low- junction emitter 

(OCI-HLE) solar cell. Eor the particular structure studied, under AMI or 

AMO illumination, this suppression increases considerably over the 

+ 

values seen in a conventional n -p solar cell of the same base doping 
concentration. Larger increases in are to be expected. These expecta- 
tions are discussed together with the advantages inherent in this type of 
HLE solar cell. 


PkQZ 8L#i'4K BU® 
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THE STRUCTURE: As is shown In Fig. 1» the H-L junction in the emitter 

results from an electron accumulation layer induced in the n-type epitaxial 

emitter region by a positive oxide charge Q^. The charge Is created near 

the Si-SiO^ interface during a low temperature oxygen heat treatment [6-8). 

To assure a good ohmic contact, a shallow n diffusion is made under the 

metallized portion of the top surface area. This n diffusion also provides 

a high- low-barrier and a small effective surface recombination velocity [9] 

17 -3 

for holes. The p-type base has a doping concentration «AA 5x10 cm , 
which approximately minimizes J_ in a cell without a back-surface field [10]. 

D 

For good cell performance, the n-type epitaxial region must satisfy four 

conditions: (a) it must be thick enough and doped heavily enough to provide 

the sheet conductance needed to make the series resistance negligible; 

17 ”3 

(b) its doping concentration must be low enough (less than about 10 cm ) 
to avoid the effects of energy band gap narrowing and severe hole lifetime 
degradation; (c) its hole lifetime must be long enough that nearly all the 
holes injected from the base diffuse through to the surface rather than 
recombining within the region; and (d) the surface recombination velocity 
under the oxide should be lov7. 

The mechanisms underlying the suppression of are simple. As the 
injected holes diffuse into the electron accumulation layer, the partially 
reflecting boundary condition [9] at the H-L junction depresses the hole 
density so that the hole recombination current is controlled by the elec- 
tron concentration N near the surface and the hole surface recombination 


velocity S^: 


dp (qn^Sp/Ng)[exp(qV/kT) - 1] 


( 1 ) 
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Here Sp Is the surface reeoinblnation velocity at the Sl-SlO^ Interrace 
covering the non-metallized portion of the device surface area; S„ oiin he 

r 

orders of magnitude less than the surface recombination velocity at an Ohmic 
contact* As (1) suggests, if the ratio Sp/Ng is small enough, Jj^ will be 
suppressed well below J„. 

D 

The oxlde-eharge- induced H-L structure of Fig, 1 is an attractive 
realization of the HLE solar cell because it avoids high conGentratlons of 
donor impurity atoms in the thin H layer near the surface. This eliminates 
heavy-doping effects [3] such as band-gap narrowing and hole lifetime 
degradation in the H layer and it may help the achievement of low values of 
Sp. The avoidance of a diffusion step to form a highly-dpped H layer 
eliminates, over the non-netalllzed portion of the surface area, the diffusion 


of unwanted Impurities and vacancies into the silicon along with the donor 


atoms. From the viewpoint of cell performance, these properties suggest 
that the OCI-HLE solar cell can exhibit good blue response in addition to 


Increased V, 


OC* 


CELL FABRICATION: Several runs of test cells were made that demonstrate 


the desired suppression. In one such run, the starting material 

£i 


consisted of a lOym thick n-type layer of resistivity =0.1 fiem grown 


epitaxially on a p-type substrate of 


thickness and resistivity 


= Q.l fiem. 


>ase 


The wafer was oxidized in dry ©2 with 0.3% TCE at 1100*C to grow a 2500A 

O 

thick oxide layer, which was later etched to 110QA to improve the anti- 
reflection properties. The temperature of 11O0'*G was chosen for this test cell 


to assure a good -quality oxide and a low value of Sp. After oxidation the 
wafer was cooled from 110O“C to 10OO*G at the rate of SG^C per hour and then 
cooled to 700 C in three hours. Holes were opened In the oxide for an n 
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contact diffusion done at 900**C for 20 minutes. The device was then heated 
In dry oi^gen for 2 hours at 700"C to increase the oxide charge density 
[8-81. 

A second run employed similar processing but used starting material with 

p . ~ 1.4 12cm and pL ~ 0.14 Hem. 
epl Hiase 


EKPBRIMENTAL DEH0NST8ATI0N OF SUPPRESSION AND LARGE Since the oxide 

E UL 

charge Q Is near the SI-SIOa interface, the total charge accumulated in 
o z 

the silicon is Qe-Qn* electron distribution is described by MOS theory 

S Q 

18 ^3 

[11]. The electron surface concentration is about 3x10 cm - 0.1 cm) 

18 -3 

and about 2x10 cm (p . - 1.4 cm). Most of the electrons in the 

6pl 

0 

accufflulatl sn layer reside within about three Debye lengths (about lOOA) from 
the surfacit. 


The ba^r diffusion length L , measured by X-ray excitation [12], was 

n 

L - 20um f I, ; devices from the first run and L 30um for devices from 
n n 

the second run. Measurement of the dark I-V characteristles showed for 
both runs that the J %a;. saturation current (determined by subtracting the 


space-charge-region recombination current (2}) was equal, within measure- 


ment error, to the base recombination current Jg calculated from knowledge 

of L . Tliis demonstrates that J- « J_. Thus is determined mainly by 
n jEi D uL 

J„ in these devices. The results of the measurements show further that the 

D 

“14 2 

emitter saturation current is suppressed to less than about 5x10 kloxsT ^ 

2 

which corresponds to a voltage, (kT/q)ln(Jg^/Jg^) > 705 mV, for Jg^ = 35 raA/cm 
(measured in these devices). Thus Is suppressed to values so low that 


it presents no barrier to the achievement of 


Vqc - 700 mV, which is the maximum 


vciliie of Vgg predicted by classical theory for a 0.1 Hem base resistivity for 


AMO 1] tuminntlon. 
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The performance of Che test cells was measured at NASA Lewis Rescarcli 
Center. The cell of the first run showed the largest V V = 634 mV 

OLf Qi> 

at 25“C. Because this test cell had much larger metal coverage on the 

2 

illuminated surface (total area = = 1.34 cm ; non-metallized area “ “ 

0.972 cm^) than the 5% metallized area of a cell optimized for AMO 
illumination, the measured should be corrected to account for the loading 
effect of the parasitic p^n junction diode under the metallized portion of 
the surface. The correction, estimated from AV^^, = (kT/q) [Jln(A.j,/Ajj) ] , is 


about 8 mV. Thus this OGI-HLE cell has an open circuit voltage of 642 mV 
at 25 “G for AMO illumination. This compares with the largest values of 


V » 610 mV previously reported for n-on^p silicon p-n junction solar 

(jlj 

cells [2,13]. 


PROJECTIONS; As was emphasized, because is completely suppressed, V 

of these OCl-HLE solar cells depends mainly on J_, which is approximately 

Inversely proportional to Cells of about 0.1 S2cm base resistivity have 

shown Lj^ of approximately 100 pro [2,14] which is much larger than the 

values of L measured for our OGI-HLE test cells. These low values of L 
n n 

may result from several origins. First, in the test cells, we grew SIO 2 
at 1100“G to assure low values of and complete suppression of at the 
possible expense of L^, Second, the starting epitaxial material had 
poor values of L^. For example, the value of L^ measured on thif epitaxial 
p-n junction prior to oxidation and other processing steps was only about 
30 pm for the material used to fabricate the test ceils of the first run. 

We are trying to increase L^ and thus Vq^ by several approaches, including 
the growth of Si02 at lower temperatures, the use of better qu-iMf/ ut .-trf; I nv. 
material, and the use of a back-surface- field [10] in conjunction with 


higher values of base resistivity. An addlvtonal improvement in can be 
expected From Improvements in the short-circuit current density which 

will be discussed below. Vdien all of the potentln] improvemeiUs .-ire 
combined, the projected maximum value of Is about 700 mV at about 25**C. 

The measured value of the fill factor FF at 25°C was 0.802. The value 
of FF increases as increases, and the projected maximum FF, corresponding 
to Vqj, " 700 mV, la about 0.83 [15]. 

2 

The best measured value of of these test cells was 35 mA/cm . 

This value resulted even though no special attention was given to achieving 
good anti-reflection (AR) properties and even though the small implies 
poor collection from the longer wavelengths of the solar spectrum. Good 
blue response was measured at a wavelength of 0.5 ym. Improvements In Jg^ 
will result from Improvements In L and the AR properties. The maximum 
is expected to occur in OCI-HLE cells made with higher base resistivity 

(about 10£2cm) and a back-surface field, as in the Hellos cell for which 

2 2 
Jg^ ■ 40 mA/cm [16]. The projected maximum is about 40 mA/cm . 

These projections correspond to an AMO maximum power conversion 
efficiency of OCI-HLE solar cells in the 17 to 18 per cent range, which 
contrasts with the maximum efficiency now seen in the silicon p-n junction 
solar cells of about 15% [16] . Because of spectral differences, the AMI 
efficiency is hi^er than the AMO efficiency. The projected maximum 
efficiency of the GCI-HLE solar cell for AMI illumination is about 20%. 

High efflcienGies also can be expected in OCI-HLE solar cells designed 


for concentrated sunlight illumination in the range of 25 to 100 suns [17]. 
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Fig, 1 Cross section of an oxide-charge- induced high-low emitter 
(OCI— HLE) solar cell. Electron accumulation layer is 
induced hy a positive oxide charge Q , 
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CHAPTER 7 DISCUSSION 

This report contains the results of research in transition. 
The first two parts of the report deal with theory and experiment 
aimfld toward an investigation of the basic mechanisms responsible 
for low vqq in n-on-p solar cells. The third part describes a new 
device structure that suppresses the emitter current and raises 
well above values seen before. Our NASA LeRC research now is 
fully focussed on the work described in the third part. In sub- 
sequent reports, we will describe in more detail the theory under- 
lying the operation of HLE devices, both OCI and diffused, and 
will give the results supplied by our ongoing experiments in HLE 


solar cells. 


